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ABSTRACT 


This article contains the new orbital elements and the 
table of the general perturbations of the X th satellite of 
Jupiter. 

A modification of Hansen’s theory suggested by Musen 
was the theoretical foundation for the programming and for 
the expansion of Hansen’s coordinates of the satellite into 
trigonometric series in four basic arguments. The com- 
plete collection of formulas is given in the exposition. Also 
described are the program and the operations for handling 
the expansions and for integrating the differential equations 
of the theory by means of iteration. This program was de- 
veloped and the actual numerical computations were per- 
formed by Charnow and Maury. The new set of elements 
represents the observations in the interval 1938-1967 better 
than the previous sets. The residuals are now of the order 
of only a few seconds of arc. Several second order effects— 
for example, the planetary perturbations and the effect of 
the variability of the orbital elements of Jupiter— will be 
treated by the subsequent work. 
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APPLICATION OF HANSEN’S METHOD 
TO THE X ,h SATELLITE OF JUPITER 


by 

Peter Musen and Jesse L. Maury 
Goddard Space Flight Center 
and 

Milton Charnow 

Computer Sciences Corporation 


INTRODUCTION 

We have accumulated the observations of the X th satellite of Jupiter published over an inter- 
val of 29 years. With this observational material, an attempt can now be made to develop a better 
theory and a better system of mean orbital elements. 

The X th satellite of Jupiter was discovered in 1938 by Seth Nicholson at the Mount Wilson 
Observatory. It belongs to the group of outer satellites of Jupiter whose motion is affected pre- 
dominantly by solar perturbations but also, to some extent, by the disturbing action of Saturn. The 
satellite is very weak visually: of 19 th stellar magnitude. Its orbital inclination toward the or- 
bital plane of Jupiter is large, about 27°. Thus the development of an analytical theory for this 
satellite is not an easy problem. Soon after the discovery, two preliminary orbits were computed; 
one by Wilson (1939) without taking into account the solar perturbations, and one by Herget (UAI) 
Circ. 727) using a modified Laplacian method. Herget included the solar effect in the computa- 
tions from the very start, and thus his preliminary orbit is a better one. Herget (1947) also sug- 
gested a modification of Encke f s device to compute the perturbed rectangular coordinates of the 
satellite. Using this method Herget, up to 1943, and later Musen produced the geocentric ephem- 
erides of the X th satellite. 

A preliminary orbit correction obtained by Musen at Cincinnati was used by Herget as a foun- 
dation for the numerical integration of the differential equations of motion of this satellite. Within 
two years after the completion of the numerical integration, Musen found that his corrected ele- 
ments were not accurate enough. This might be explained by the mutual influence between the 
correction of elements and the perturbative effects. That influence could not be investigated then, 
because of lack of technical means. The possibility of such investigation exists now and enables 
us to obtain better elements as well as better general perturbations. 

Lemekhova (1961) has applied the theory of Delaunay (1860, 1867) to the motion of Jupiter X 
and obtained a corrected set of elements using observations from 1938 to 1942. She reduced these 
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observations to 14 normal places. The largest residuals after her orbit correction still remain 
36” in right ascension and 24” in declination. 

The theory of Delaunay was used by Lemekhova because it was the only analytical lunar theory 
extensively developed. One must bear in mind, however, that Delaunay had in view the application 
of his theory to the moon and not to the satellites of Jupiter. For this reason, he did not include 
the higher powers of the inclination in his expansions. It is questionable whether his theory can 
be applied uncritically to a satellite with the high orbital inclination of Jupiter X. In order for the 
theory of Delaunay to be applied to Jupiter X too, it must be extended several orders higher. This 
is a formidable task that has been undertaken by several groups, but not yet completed. So far no 
extension of the theory of Delaunay exists that would enable us to obtain an accurate representa- 
tion of observations of such satellites as Jupiter X. The authors, however, have no doubt that in 
the future such an extension will become a reality. 

The work of Lemekhova still retains much of its importance because it permits us to form an 
idea of the relative importance of terms with different arguments. During Musen’s stay in Cin- 
cinnati, Herget called his attention to the fact that the lunar theory of Hansen (1862) can be applied 
to obtain the expansion of the solar perturbations of the outer satellites of Jupiter. Since then 
several possible modifications of Hansen’s theory have been developed by Musen (1963, 1965) at 
Goddard Space Flight Center in a form adaptable to the use of electronic machines. One of those 
modifications (1963) is being utilized in this work. The programming and the numerical computa- 
tions were performed by M. Charnow (1966) and Jesse Maury. 

Hansen’s lunar theory can be applied to the cases of high inclination. It is also usable for 
obtaining the expansions of perturbations up to any degree of accuracy compatible with the ac- 
curacy of observations. These expansions are obtained in the form of trigonometric series with 
purely numerical coefficients, and this makes the programming easier. So far this is the only 
theory which permits relatively rapid calculations. Some numerical inaccuracies do exist in 
Hansen T s theory of the moon. However, his basic theoretical idea is correct and can be applied to 
satellites on which solar effects strongly dominate and for which several small non- gravitational 
effects, important in the lunar theory, can be neglected in the computation of yearly ephemerides. 
The intermediary orbit of Hansen T s theory is a rotating ellipse of fixed shape. The use of the var- 
iational solution (Hill, 1878) is essential only when a literal development in powers of the constants 
of the eccentricities and of the inclination is being pursued (Brown, 1897-1908). The variational 
part of the theory for the satellites is not necessarily the most significant part numerically. The 
largest amplitude normally is associated with the evection and contains the eccentricity of the sat- 
ellite orbit as a factor. Thus, from the numerical point of view, the use of the mean rotating 
ellipse as an intermediary represents a natural start. 

The appearance of squares of small divisors in the process of double integration presents 
considerable difficulty in every lunar theory. In Hansen’s theory, the squares of small divisors 
are associated with the long period terms in the perturbations of the mean anomaly and with the 
elements. The amplitudes of these terms are small. In the process of iteration they are obtained 
as a ratio of two small numbers. Thus a loss of accuracy is inevitable. Every lunar theory has a 
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special procedure for handling these terms. Here we made use of Hansen T s transformation (1862) 
of the element E. The expression for this element is written by Hansen in such a manner that the 
terms to be integrated twice contain small factors. The presence of these factors diminishes the 
influence of division by the mean motion of long period arguments in the double integration. 

After the programming had been completed, we applied several checks to test its correctness. 
The expansion of the disturbing function pertaining to the first cycle of integration must coincide 
with that obtained by Delaunay when the numerical values of the eccentricities and of the inclina- 
tion (both constants) are substituted. We found that the agreement was very good and was always 
in accord with the accepted numerical accuracy. The process of iteration was checked by the 
partial reproductions of Hansen’s lunar theory. The final output of the theory comprises the 
Fourier expansions of the Hansen coordinates and of the Euler parameters of the satellite. Their 
accuracy as given here is quite sufficient to produce a good ephemeris of the X th satellite. 

Using Lemekhova’s elements and our expansions of perturbations, M. Charnow and J. Maury 
computed an ephemeris for 1967. The satellite was found by E. Roemer in very close proximity 
to the predicted position. Lemekhova in her work expressed the opinion that some observations 
of Jupiter X still remain unpublished. The authors also hold this not impossible. It would be highly 
desirable that such unpublished observations be made available to the scientific community. 

The orbit correction and the perturbations have now been recomputed using the new observa- 
tions. The present work contains the final results of these computations, the complete set of 
formulas used, and also the comparison with the observations. One can see that the set of new 
elements represents the observations much better than the set of elements with which we started. 

We prefer the general perturbations because they are superior to numerical integration in the 
description of physical characteristics of an oscillatory system. All the theoretical and technical 
difficulties associated with the work on the general theories are compensated for by the gratifying 
results. We obtain an accurate representation of the observations, an accurate prediction, and 
also an insight into the behavior of the orbit over a long interval of time. 

Work on the general theories of motion of bodies in the solar system must proceed because 
flights to the planets and to their satellites will become a reality in the future. Some day we will 
have artificial planetary satellites too. The numerical and theoretical techniques must be ready 
to meet these new challenges. 


COLLECTION OF FORMULAS 

The notations used in this work are the same as in Musen’s previous theory (1963). 

A Bessel-function routine is employed to generate the numerical coefficients of the input 
series (p/ a 0 ) cos 0 through a'/r'. These series have literal arguments. The major iterative loop 
begins with the calculation of the s. series. The disturbing function is generated from the s and 
the p series. When the series representing the d^/dt, d\ 2 /dt, and d\ 3 /cit functions have been 
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formed, the values of n 0 y, n 0 a, and n 0 v are determined by imposing the condition that the series 
representing those functions must contain no constant terms. 

The values n 0 y, n 0 a, n 0 77, and n 0 y' relate to the motion of the lunar perigee, lunar node, 
and solar perigee. The series 

h ol 

IT ’ [T] , [¥] , [A..] . (i = 1 , 2, 3 , 4 ) 

are obtained by formal integration. The n 0 Sz series is computed from the dn 0 Sz/dt series after 
c 1 and c 2 have been calculated, so that n Q Sz contains neither constant nor sin g terms. 

The n 0 Sz series has been computed; now the 

^0 h _ 

"h ’ B , ^ ^ . w and v 

are constructed for final output or for use in the next iteration. They are complete when the values 
of n 0 y , n 0 a, and n Q 77 have converged. 
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Designating by 


T? ’ [T] , M , [\J (i = 1 , 2 , 3 , 4 ) 


the series obtained by formal integration, we have 
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DETERMINATION OF THE LONG-PERIOD EFFECTS 


The terms containing the short period arguments g or g' can be handled by the process of 
iteration easily. The difficulty arises, for reasons explained in the Introduction, when we want to 
determine the long period terms in the perturbations of the mean anomaly, among which the term 
with the argument 2<u is the most important. For the computation of these terms, a special pro- 
cedure is necessary. The most significant long period terms in the perturbations of the mean 
anomaly are transferred into n Q Sz through the element H. An additional integration is performed. 
Consequently, we need a form of S split into two parts. The first part will be affected by a single 
integration only. The second part, which will be affected by a double integration, will contain 
small factors which will diminish the effect of division by the mean motion of the long period ar- 
guments. We are here using, without modification, Hansen's formulas given on page 374 of the 
Darlegung (Hansen, 1862): 


3a o n " 3 TT 0 ~ e o 2 ' V 1- e ' 2 (lT cos 1 + k ) (fr) 




^ -l) + f (l-e 0 ’)(T» + *) +Z 


cos I = 1 - 2\j 2 - 2V 2 2 , 



77 ^- = - 3 ( a o n i) - 4 ( a o n 2 ) - s (a 0 n 3 ) 


The series 


" 1 = - a + a 2 - a 2 + • • • - 1)2 , t 2 . r 

can be taken from the previous iteration. 

At each iteration step we shall use these formulas only for the purpose of determining the 
long period terms. In this case we can disregard the constant of integration k. 
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COMPUTATION OF THE COORDINATES AND VELOCITIES 


After the process of iteration has been completed, we determine the components of the satel- 
lite position and velocity vectors from the following system of formulas: 


Aj (a) 
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0 0 
+ cos a - sin a 
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Let the Gibbsian vectors of Jupiter's orbit be 


P\ Q\ S' . 


The components of these vectors and of the satellite position and velocity vectors are referred to 
the mean equator and equinox 1950.0. 

We have 
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The geocentric right ascension and declination and the distance from the earth are computed from 
the standard formulas 


p cos 8 cos a = x + x' + X , 

p cos 8 sin a = y + y ' + Y , 

p sin 8 = z + z ' + Z , 

where x', y', z' are the heliocentric coordinates of Jupiter and X, Y, z the geocentric coordinates 

of the sun. 

ORBIT CORRECTION 

The formulas used for the orbit correction represent a transformation of the Eckert-Brouwer 
method (1937). In this transformation we make use of our expressions for the position and velocity 
vectors in the inertial system. 
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PROGRAMMING 

The adaptation of Hansen's lunar theory as a model has made possible the effective applica- 
tion of digital computing equipment to this problem. Hansen required a qualitative personal judg- 
ment at every step of the development of his long Fourier series. We replace this with a feedback 
controlled iteration which now can be automated. A computer program can be constructed of a 
sequence of relatively short steps: this sequence is repeated until we reach the accepted numerical 
accuracy. These short steps are comparatively straightforward to program and easily checked. 

The basic tools of the Hansen lunar theory are Fourier series with five periodic arguments, 
g, g’, o, a>‘ , and y. Since all needed operations on these series (series multiplication, series ad- 
dition, differentiation, integration, etc.) are closed, a set of computer subroutines to execute these 
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operations can readily be defined. In the technique adopted by the authors, each term of a series 
is composed of two words. The first is the numerical coefficient in the usual floating-point form. 
The second is a "logic” word containing the periodic arguments. This logic word relies on posi- 
tional representation and bias (normalization) around 50. It consists of a base logic word having 
the value +005050505050505050 with the integral multipliers of the periodic arguments introduced 
into it in the order g, g f , co 9 a>‘ , y . The coefficient of the argument g can assume any integral value 
between 0 and 99. The coefficients of the four remaining arguments can assume any integral value 
between -49 and +49 normalized by 50. By maintaining all logic words so that the first non-zero 
argument is positive (i.e., A sin (-6) = - A sin (0) and B cos (-0) = B cos (0) ), it was possible to use 
the sign of the logic word to denote sine or cosine. This effects a saving in computation time and 
storage space. Cosine is denoted by a positive logic word; sine by a negative word. An example 
is +029329805 sin (lg - 2g # + 2co - 2^'), the evection term in n 0 Sz, In the computer, this term is 
represented by the numeric, floating-point coefficient +029329805 followed by the logic word 
-014852485050505050. 

The constant term in any series is carried as a coefficient of cos (0): that is, a floating-point 
coefficient followed by the base logic word +005050505050505050. This representation is unique. 

To each series is prefixed a single word denoting the number of terms in the series. The terms 
of a series are arranged in descending order of the absolute value of the numeric coefficients. 

This representation of series was used to program a complete set of computer subroutines 
for the algebraic manipulation of Fourier series. This package of subroutines can multiply, add, 
subtract, scalar multiply, differentiate, integrate, and evaluate trigonometric series, as well as 
perform several specialized functions such as extracting a specific term with a given argument. 
Series multiplication employs the standard half-angle formulas: addition and subtraction of the 
arguments are realized through logical addition and subtraction of the logic words (Charnow, 1961): 
all five arguments are handled simultaneously without recourse to breaking apart the logic word. 

In executing any of these operations, all duplication of terms is eliminated (maintaining a positive 
non-zero first argument is an aid in this). An operation on a series or a binary operation on two 
series is terminated by the exhaustion of terms; or, since the terms of a series appear in descend- 
ing order, it may be terminated on the basis of a numerical criterion applied to the coefficients 
of the resultant series. At the conclusion of any operation which may produce a disordered series, 
the resultant series is ordered by a separate subroutine, the arranger. By the use of these tools, 
algebraic series operations are programmed almost as easily as numerical calculations. 

The original version of this program was written in a pseudo-machine language (Maury, 1964), 
(Gorman, 1964) on the IBM 7090 and 7094. This version was capable of operating on series with 
fifty or fewer terms. Subsequently, versions allowing longer and longer series were developed. 

The goal was to ensure that the computation would be controlled only by numerical criteria and not 
by storage restrictions. To realize this, the program was shifted to the Univac 1107 to take ad- 
vantage of the high-speed, random-access drums. These drums were used as temporary storage 
for series, while the bulk of the computer storage was reserved for two operand registers of 
10,000 words each and a result register of 20,000 words. 
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As the length of the series increased, the time to perform the series operations also increased 
(exponentially in the case of the series multiply operation). To ease the time burden, a number of 
the subroutines were rewritten in machine language and the program was again shifted, this time 
to the Univac 1108. Currently, all the main subroutines including multiplication, addition, sub- 
traction, and the arranger are in assembly language, while the statements reflecting the lunar 
theory are in symbolic language. This version of the program can handle series having up to 
2,500 terms, and 750,000 words of high-speed (4.25 ms access time) drum storage is reserved for 
storing series. The computer storage contains the subroutines, the lunar theory, a small resident 
I/O monitor, and the 40,000 words comprising the three series registers. 

The application of this program to the X th satellite of Jupiter required two hours. The first 
iteration took only 40 seconds, since the perturbations are zero during this iteration. The second 
iteration required 16 minutes. The third iteration, in which all the perturbations were completely 
represented and the series had just begun to converge, required 26 minutes. Thereafter, the time 
to complete one iteration dropped off markedly, to approximately 13 minutes. Convergence was 
achieved in eight iterations. 


CORRECTED ELEMENTS OF THE X th SATELLITE 

The initial elements of the motion of Jupiter X are derived from those developed by E. N. 
Lemekhova: 

Epoch: 1938 July 27.3128 = J.D. 2429106.8128 


X Satellite Jupiter 


The elements are referred to the 
orbital plane of Jupiter 

Mean ecliptic and equinox 1950.0 

g 0 = 216.°6928 

K‘ 

= 329°20254 

n = 1 “384557 per day 

n ' 

= 0°083091 per day 

a = 0.078345 a.u. 

a' 

= 5.202561 a.u. 

e = 0.10739 

e ' 

= 0.048398 

= 244^7381 

co' 

= 273°57223 

cj 0 ' = 113.5835 

n' 

= 99.92939 

I = 27.5748 

i ' 

= 1.30614 


The system was adjusted for use with a modified Hansen instead of Delaunay theory. 
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Using the initial elements, the Hansen parameters n 0 y, n 0 a, and n 0 77 were obtained from the 
general perturbation theory: 

n 0 y = 1?64886 per year 
n 0 a = 1.15605 per year 
n Q V = 0.06700 per year 

The mean motions of the perigee and node are : 

mean motion of perigee = n 0 (y - 2-n) = +lf51486 per year 
mean motion of node = - n 0 (a. + 77 ) = -l.°22305 per year 

We obtained the following corrections: 

Ag 0 = +3 f 6449 ± 0?0310 
A x = -3.4020 ± 0.0308 
sin IA co 1 — - 0.1320 ± 0.0090 
Al = +0.0256 ± 0.0072 
Ae = +0.00851 ± 0.00009 
Aa/ a = -0.000128 ± 0.000001 


The system of the corrected elements is 

Epoch J.D. 2429106.8128 
g 0 = 220^3377 
n = l.°384687 per day 
a = 0.078335 
e = 0.11590 

co 0 = 241?0835 

« 0 ' = 113f2985 

I = 27.°6004 
n 0 y = 1^65492 per year 

n 0 a = lfl 6939 per year 

n Q 77 = 0.°06831 per year 

This set of elements was obtained through the repeated use of a differential correction pro- 
cedure. The new elements resulted in better agreement between the observed and predicted posi- 
tions of the satellite. 


s 
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The Representation of Observations Used. 

Observed O - C 





a (1950.0) 

8(1950.0) 

Aa cos 8 

AS 


1 

1938 July 

6^3743 

22 h 17 m 15 s .21 

-11°13 , 32V3 

-173 

+276 

Mt. Wilson 

2 


6.4667 

22 17 13.69 

-11 13 41.6 

-0.9 

+1.5 

Mt. Wilson 

3 


9.4229 

22 16 22.89 

-11 18 23.4 

-1.5 

+1.4 

Mt. Wilson 

4 


9.4521 

22 16 22.38 

-11 18 26.0 

-0.9 

+1.8 

Mt. Wilson 

5 


27.3361 

22 09 13.67 

-12 00 14.9 

-0.2 

+0.0 

Mt. Wilson 

7 


28.3174 

22 08 45.05 

-12 03 06.3 

-0.3 

+0.3 

Mt. Wilson 

8 


29.4667 

22 08 10.91 

-12 06 30.2 

+0.1 

+1.1 

Mt. Wilson 

10 

Aug. 

24.4167 

21 54 07.40 

-13 33 16.2 

+2.4 

+0.2 

Mt. Wilson 

11 


25.2129 

21 53 41.88 

-13 36 00.3 

+3.2 

-0.8 

Mt. Wilson 

12 

Oct; 

18.1285 

21 39 13.15 

-15 25 53.6 

+0.7 

+0.1 

Mt. Wilson 

13 


20.1257 

21 39 25.45 

-15 25 54.0 

+0.7 

-0.5 

Mt. Wilson 

14 


23.1319 

21 39 50.17 

-15 25 15.1 

+0.2 

-0.0 

Mt. Wilson 

15 

Nov. 

21.1188 

21 49 30.06 

-14 41 12.0 

+1.0 

+0.2 

Mt. Wilson 

16 

1939 July 

15.4205 

0 31 26.11 

+ 1 36 35.3 

+3.0 

-0.6 

Mt. Wilson 

17 


15.4528 

0 31 26.49 

+ 1 36 36.4 

+2.2 

-0.7 

Mt. Wilson 

18 

Aug. 

16.4365 

0 33 47.03 

+ 1 33 14.1 

+1.3 

+2.6 

Mt. Wilson 

19 


16.4806 

0 33 46.72 

+ 1 33 11.8 

+1.0 

+2.4 

Mt. Wilson 

20 

Oct. 

8.3998 

0 15 53.39 

+ 0 00 48.7 

-0.3 

+3.3 

Mt. Wilson 

21 

Dec. 

15.0885 

0 01 26.93 

- 0 48 17.8 

-2.8 

-0.6 

Mt. Wilson 

24 

1940 Sept. 

8.3736 

2 57 26.16 

+16 02 05.9 

-0.2 

-0.1 

Mt. Wilson 

25 


8.4438 

2 57 25.74 

+16 02 03.0 

+0.1 

-0.6 

Mt. Wilson 

26 

Oct. 

25.3488 

2 38 37.74 

+14 14 20.0 

+0.9 

-0.7 

Mt. Wilson 

27 


25.3731 

2 38 36.88 

+14 14 15.2 

+1.6 

-0.3 

Mt. Wilson 

28 

Nov. 

1.2160 

2 34 23.98 

+13 49 17.5 

+0.0 

-1.0 

Mt. Wilson 

29 


1.4181 

2 34 16.35 

+13 48 33.3 

+0.1 

-0.6 

Mt. Wilson 

30 

1941 Dec. 

23.2569 

4 56 06.78 

+22 24 12.2 

-1.0 

+0.2 

Mt. Wilson 

31 

1942 Feb. 

17.2342 

4 44 28.41 

+22 23 09.5 

+0.7 

-1.1 

Mt. Wilson 

32 


18.2618 

4 44 36.92 

+22 23 19.4 

+2.0 

-1.5 

Mt. Wilson 

33 

Nov. 

8.4951 

7 53 12.49 

+21 30 26.6 

-0.1 

-0.4 

Mt. Wilson 

34 


9.5090 

7 53 17.86 

+21 30 12.4 

+0.2 

-0.6 

Mt. Wilson 

35 

1943 Jan. 

6.2400 

7 32 51.54 

+21 57 36.7 

+6.1 

+1.1 

Mt. Wilson 

36 


6.3230 

7 32 48.03 

+21 57 40.5 

+5.6 

+0.5 

Mt. Wilson 

37 

1951 Sept. 

30.29101 

0 39 41.51 

+ 2 17 48.3 

-1.1 

-3.8 

Mt. Wilson 

38 


30.30663 

0 39 41.08 

+ 2 17 47.1 

-1.6 

-3.3 

Mt. Wilson 

50 

1958 Apr. 

25.30100 

13 34 27.90 

- 7 41 04.4 

-2.2 

-2.0 

Flagstaff 

51 


25.36887 

13 34 26.09 

- 7 40 53.6 

-1.3 

-2.3 

Flagstaff 

52 

1967 Jan. 

7.31247 

8 17 49.28 

+20 23 30.2 

-0.1 

+4.2 

Tucson 

53 


7.38122 

8 17 47.19 

+20 23 34.4 

+1.4 

+3.3 

Tucson 

54 

Feb. 

11.17166 

7 57 31.87 

+21 06 02.6 

-6.3 

+2.7 

Tucson 

55 


11.25846 

7 57 28.83 

+21 06 08.1 

-7.3 

+2.9 

Tucson 
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+ 0.00094 

1 

-3 

+1 

-3 

- 0°007 

+ 0.00005 

2 

-4 

+3 

-3 

- 0°001 


0 

+1 

-3 

+1 

+ 0.007 


1 

-3 

+2 

-2 

+ 0.213 

- 0.00170 

2 

-4 

+4 

-4 

+ 0.015 

- 0.00010 

0 

+1 
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+0 

+ 0.043 

- 0.00002 
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+3 

-3 
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-3 
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-2 
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0 
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+2 
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-2 

- 0.003 
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0 

+1 

+0 

+2 
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-1 
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-2 
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-2 
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0 
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+1 
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1 

-2 
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-2 
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-2 
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+0 
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- 0.002 
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-1 
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- 0.070 

+ 0.00058 

2 

-1 

+2 

+0 


- 0.00001 i 

0 

+2 

-2 

+2 

+ 0.040 

- 0.00074 

1 

-1 

+2 

-2 

- 0.060 

+ 0.00051 

2 

-1 

+3 

-1 

+ 0.001 

j 

0 

+2 

-1 

+1 

+ 0.008 


1 

-1 

+2 

+0 

+ 0.011 

- 0.00009 

2 

+0 

+0 

+0 


+ 0.00003 j 

0 

+2 

+0 

+0 

- 0.006 


1 

-1 

+3 

-1 

- 0.001 

+ 0.00001 

2 

+0 

+2 

+0 

- 0.027 

+ 0.00060 ! 

0 

+2 

+0 

+2 

- 0.476 

+ 0.00037 

1 

-1 

+4 

-2 

- 0.001 

+ 0.00001 

2 

+0 

+4 

+0 

+ 0.019 

- 0.00013 

0 

+2 

+1 

+1 

- 0.002 


1 

+0 

-2 

+0 

- 0.022 

+ 0.00019 i 

2 

+1 

+0 

+0 

- 0.001 

+ 0.00003 . 

0 

+2 

+2 

+0 

+ 0.001 


1 

+0 

+0 

+0 


+ 0.00039 1 

2 

+1 

+1 

+1 


-o.ooooi : 

0 

+2 

+2 

+2 

+ 0.007 


1 

+0 

+2 

-2 

- 0.004 

+ 0.00003 

2 

+1 

+2 

+0 

+ 0.002 

- 0.00002 

0 

+3 

-2 

+2 

+ 0.009 

- 0.00009 

1 

+0 

+2 

+0 

- 1.861 

+ 0.01628 

2 

+2 

+0 

+2 

- 0.001 

+ 0.00002 | 

0 

+3 

-1 

+3 

- 0.004 

+ 0.00001 

1 

+0 

+4 

+0 

+ 0.004 

- 0.00004 

2 

+2 

+2 

+2 

+ 0.002 

- 0.00002 

0 

+3 

+0 

+2 

- 0.054 

+ 0.00006 

1 

+1 

-1 

+1 

+ 0.001 

- 0.00001 

3 

-5 

+4 

-4 

+ 0.002 

- 0.00002 

0 

+3 

+1 

+3 

- 0.001 


1 

+1 

+0 

+0 

- 0.058 

+ 0.00054 

3 

-4 

+4 

-4 

+ 0.008 

- 0.00008 

0 

+3 

+2 

+2 

+ 0.001 


1 

+1 

+0 

+2 

+ 0.002 

- 0.00002 

3 

-3 

+2 

-2 

+ 0.004 

- 0.00008 

0 

+4 

-2 

+2 

+ 0.003 

- 0.00001 

1 

+1 

+1 

+1 

+ 0.024 

- 0.00022 

j 3 

-3 

+3 

-3 


- 0.00001 . 

0 

+4 

-2 

+4 

+ 0.002 

- 0.00001 

1 

+1 

+2 

+0 

- 0.012 

+ 0.00011 

i 3 

-3 

+4 

-2 

- 0.001 

+ 0.00001 

0 

+4 

-1 

+3 

+ 0.001 


! l 

+2 

-2 

+2 

+ 0.001 

- 0.00001 

3 

-2 

+2 

-2 

+ 0.021 

- 0.00044 

0 

+4 

+0 

+2 

- 0.004 


: 1 

+2 

+0 

+0 

- 0.002 

+ 0.00002 

3 

-2 

+4 

-2 

- 0.010 

+ 0.00009 

0 

+4 

+0 

+4 

- 0.001 


; l 

+2 

+0 

+2 

- 0.052 

+ 0.00048 

; 3 

-1 

+2 

-2 

- 0.001 

+ 0.00001 | 

0 

+5 

-2 

+4 

+ 0.001 


l 

+2 

+1 

+1 

+ 0.001 

- 0.00001 

! 3 

+0 

+2 

+0 

- 0.001 

+ 0.00005 1 

1 

-5 

+2 

-2 
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l 
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+2 
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1 
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CONCLUSION 


It has been shown that electronic data processing machines can be used to automatize the pro- 
duction of semi- analytical theories of the "natural” satellites. The theory of motion of the X th 
satellite of Jupiter as given in this article represents the solution of the so-called "main problem." 
The direct solar effects were expanded under the assumption that the heliocentric motion of Jupiter 
is Keplerian. We are justified in such an approach because the solar perturbations are undoubtedly 
strongly dominant in the motion of Jupiter X. This can be seen from the fact that the new set of 
elements leads to a considerably better representation than the old one. The difference between 
the observed and the computed positions is now reduced to only a few seconds of arc. However, to 
make the theory more complete, several other secondary small effects must also be considered. 

We have still to include in the theory the direct planetary influence as well as the influence of the 
variability of the elements of Jupiter on the motion of the satellite. We expect these effects to 
change the difference between the observed and the computed position slightly, about 1" - 2". 

The comparison between the theory of Hansen and that of Delaunay also is a topic of consider- 
able theoretical and numerical importance. The authors plan to expand and to discuss these prob- 
lems in their next paper. 
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